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We study electron-phonon interaction induced decoherence between two-electron singlet and 
triplet states in a semiconductor double quantum dot using a spin-boson model. We investigate 
the onset and time evolution of this dephasing, and study its dependence on quantum dot param- 
eters such as dot size and double dot separations, as well as the host materials (GaAs and Si). At 
the short time limit, electron-phonon interaction only causes an incomplete initial Gaussian decay 
of the off-diagonal density matrix element in the singlet-triplet Hilbert space, a complete long-time 
exponential decay due to phonon relaxation would eventually dominate over two-spin decoherence. 
We analyze two-spin decoherence in both symmetric and biased double quantum dots, identifying 
their difference in electron-phonon coupling and the relevant consequences. 



I. INTRODUCTION 

Significant experimental progresses in the study of 
semiconductor spin qubits in the past few years 1-1 7 have 
reconfirmed the confined electron spins as one of the lead- 
ing candidates for the building block of a solid state quan- 
tum information processor. A decade of theoretical stud- 
ies has mostly clarified single spin decoherence channels 
and their relative importance in semiconductor quantum 
dot (QD) and donor confined electrons, 18-40 with hyper- 
fine interaction to the lattice nuclear spins as the main 
culprit for electron spin decoherence. 

Decoherence of two-spin states in a coupled dou- 
ble quantum dot is crucial to the operation and 
scale-up of exchange-based spin quantum computer 
architectures. 41-46 Since nuclear spins are the main 
sources of single spin decoherence in GaAs quantum dots, 
where most experimental progress has been made, exist- 
ing theoretical studies have focused on the decohering ef- 
fects of the nuclear spins. 47-50 In addition, since exchange 
coupling is electrostatic in nature, exchange-coupled elec- 
trons are vulnerable to charge noise and other orbital 
fluctuations that have an electrical signature. 48,51-55 For 
example, we have shown how gate noise 51 and back- 
ground charge fluctuations 52 lead to pure dephasing by 
introducing noise into exchange splitting of a double dot. 

Electron-phonon interaction is intrinsic to any solid 
state system, 56,57 and semiconductor nanostructures are 
no exception. It is therefore important to consider the 
role of electron-phonon interaction in electron spin de- 
coherence. While electron-phonon interaction is gener- 
ally not spin-dependent, it can affect spins when com- 
bined with other interactions. For example, in a sin- 
gle quantum dot, electron-phonon interaction can assist 
single-electron spin flip or two-spin transitions in combi- 
nation with spin-orbit interaction 18,24,27,58-63 or hyper- 
fine interaction. 20 ' 47,50 In the case of donors, the strongly 
localized electron wave function and the resulting lattice 
strain lead to a direct spin-lattice interaction, so that 
electron-phonon interaction can cause pure dephasing for 
a single spin. 64 For a pair of donors close to each other, 



the two-electron spin states can be mixed by the hyper- 
fine interaction with the P nuclear spins, which allows 
two-spin relaxation via phonon emission. 65 



In this work we study the decoherence effects of 
electron-phonon interaction on two-electron-spin states 
in semiconductor double quantum dots (DQD's). Singlet 
and triplet states are two-spin eigenstates for exchange- 
coupled electrons in the absence of spin-orbit interaction 
and inhomogeneous magnetic fields (otherwise electron- 
phonon interaction can lead to relaxations between sin- 
glet and triplet states 50,62,63,65 ). These two types of 
states have different charge distributions because of their 
different spin symmetry. We show that this difference 
in electron charge density distribution leads to different 
dressing by the phonons, without the involvement of the 
excited states and/or spin-orbit interaction. This differ- 
ence in phonon dressing then leads to pure dephasing 
between singlet and triplet states. The systems we con- 
sider include coupled quantum dots in GaAs and Si, both 
regarded as promising candidates for qubits in spin-based 
quantum information processing. 



The paper is organized as follows. In Sec. II, we in- 
troduce the electron-phonon interaction in GaAs and Si. 
Combined with knowledge of two-electron states, we ob- 
tain the effective interaction Hamiltonian in the form of 
a spin-boson model, and clarify the dynamics of two-spin 
dephasing. In Sec. Ill, we present our results, quantify- 
ing the time scale of two-spin dephasing in both GaAs 
and Si, and in both symmetric and biased double dots, 
and identifying the most important types of electron- 
phonon interaction. Finally, in Sec. IV, we discuss the 
implications of our results on spin and exchange-based 
quantum information processing, and we give our con- 
clusions. 
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II. THEORETICAL FORMALISM 
A. Electron-phonon interaction in GaAs and Si 

The general electron-phonon interaction Hamiltonian 
in a semiconductor takes the form 57 

Hep = ^A(q)p(q)(a q ,A + a f _ q J , (1) 

q,A 

where a q ,A and al q A are phonon annihilation and cre- 
ation operators with lattice momentum q and branch in- 
dex A, and p(q) is the electron density operator. For 
this work we consider the electron interaction with both 
acoustic and optical phonons. 

For semiconductors with polar characteristics, such as 
GaAs and In As, electron-phonon interaction is gener- 
ally strong, including deformation potential (DP) inter- 
action and piezoelectric (PE) interaction with acoustic 
phonons, and polar (PO) interaction with longitudinal 
optical (LO) phonons. Deformation potential interaction 
in GaAs only couples electrons to longitudinal acoustic 
(LA) phonons, 



^GaAs(q)=^ 



pVojc 



(2) 



where D is the deformation constant, p is the mass den- 
sity, V is the volume of the crystal, and w q is the angular 
frequency of the phonon mode q. For GaAs D = 8.6 eV 
and p = 5.3 x 10 3 kg/m 3 . For piezoelectric interaction in 
a zinc-blende lattice, 



M, 



PE 
GaAs 



(q) = i 



pVujc 



where e is the elementary electric charge, en is an elas- 
ticity tensor component, £ is the polarization vector, and 
q is the unit vector along q. For GaAs eu = 1.38 x 10 9 
V/m. Notice that PE interaction can couple electrons to 
both LA and transverse acoustic (TA) phonons. For po- 
lar interaction with LO phonons in bulk polar materials 
such as GaAs and InAs, 



M£ a ° As (q) 



' 2ire 2 huj 



LO 



q 2 V 



1 



(4) 



where eoo and eo are the high- and low-frequency dielec- 
tric constants, and Tiwlo is the zone-center LO phonon 
energy. For GaAs, = 10.89, eo = 12.9, and Swlo = 
36.25 meV. In a quantum well with well width a z , where 
barrier materials have different dielectric constants than 
the well itself, the LO phonons are confined, so that the 
LO phonon wave vectors along the confinement direction 
can only take discrete values of q z = nn/a z , with n being 
positive integers. 66 

For Si, which has a vanishing PE interaction because 
of the inversion symmetry of its lattice, the DP inter- 
action has similar strength as in GaAs, and can couple 



electrons to both acoustic phonon branches. However, 
there is no interaction between conduction electrons and 
optical phonons in Si. 57 The conduction band of bulk Si 
has a six-fold degeneracy at its bottom, 57 so that the DP 
electron-phonon interaction takes on a more complicated 
form. 57 For an electron in a particular valley along the k 
direction, 



= E d Tr{e} + E u (k-e-k). 



(5) 



where Ed and E u are the dilation and shear deformation 
potential constants, and e is the strain tensor of the lat- 
tice due to lattice vibrations. For Si, Ed = 5.0 eV and 
E u = 8.77 eV for electrons at the bottom of the con- 
duction band. 57,67 For a two-dimensional quantum dot 
(in the xy-plane) whose electronic ground orbital state 
involves only the z and — z valleys, 



m& (q) 



i «-DP,TA/ \ 

^ S i (q) 




1 + ^ 2 ) , (6) 



(7) 



Having obtained the explicit forms of the electron- 
phonon interaction Hamiltonians in both GaAs and Si, 
we can now project them onto a specific electronic state 
basis. Below, we discuss these projections in both sym- 
metric and biased double quantum dots with two elec- 
trons. 



2eei4 {q x q y £ z + q y q z £ x + q z q x £y) , 

(3) 



B. Charge distribution of two electrons in a 
symmetric double quantum dot 

For two electrons in a DQD, the electron density op- 
erator p(q) in the general electron-phonon interaction 
Hamiltonian Eq. (1) takes the form p(q) = e* q ri + 
e iq-r 2 _56 with ^ ne knowledge of electron orbital states, 
we can calculate the matrix elements of p(q). 

When two spin qubits are exchange-coupled in an un- 
biased symmetric DQD, their orbital states are sym- 
metric or anti-symmetric if their spin state is a singlet 
(|tl - ;t>/\/2) or a triplet (\U - lt)/V2, |tt), III))- 
Within the Heitler-London approximation, the two spa- 
tial wave functions can be written as 



Ms) = 



\L{l)R{2) + R{l)L{2)) 



\/2(l + S 2 )' 

Was) = -y=l—\L(l)R(2)-R(l)L(2)), (8) 

where L and R refer to the ground single-electron orbital 
states in the two dots, S — (L\R) is the overlap integral, 
and 1 and 2 are indices for the two electrons. 

Now we can project the electron-phonon interaction 
into the singlet-triplet Hilbert space. All three triplet 
states have the same orbital wave function and cannot 
be differentiated by electron-phonon interaction. The 
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Hilbert space of interest is thus only two-dimensional, 



with the corresponding basis states 



R(1)L(2)) x -L\U - It) 
R(1)L(2)) x ^|tl + It)- 



and 



-r-± \L(\)R{2)- 

1 |£(l)i2(2) - 



V 2 (i-s 2 ) 
Since the Hamiltonian has 

no spin-dependence, the 2x2 electron-phonon interac- 
tion Hamiltonian is diagonal: 



ff eff = ^M x (q)A^(a q , A + a t _ qA ), (9) 

where cr^ is a Pauli matrix in this two-dimensional two- 
electron Hilbert space (it is not for single electron spins) , 
and the charge distribution difference A$ is given by 

Af,= \[{'<pAs\p{^)\ipAs)-{'4>s\p{ci)\ips)] =A 4> (q ll )f(q z ). 

(10) 

Here f(q z ) is determined by the z-direction (growth di- 
rection) wave function, and there is no transition be- 
tween subbands created by z-confincment. For an infi- 
nite square well with width a z and for acoustic phonons 
(whose wave vectors are not limited by the quantum well 
confinement), 



f{Qz 



sin q z a z 



q z a z (q z a z ) 2 - it 2 



(11) 



For LO phonons, q z are discrete: q z — mir/a z , with 
m = 1,2, .... In the present calculation there is no in- 
tersubband transition, so that 



f(q z = 2nir/a z ) = , 
while for q z = (2n + l)n/a z , 

'(2n+ 1)tt\ _ (-1)™ +1 



(n- l/2)(n + l/2)(n + 3/2) 



(12) 



(13) 



For a symmetric DQD, the singlet state has larger 
charge density in between the two dots, while the triplet 
has larger charge density at the far ends of the DQD. 
The resulting difference in charge distribution has a fi- 
nite electrical quadrupole moment and gives ^(qy) its 
q-dependence: 



4 ym (qn) = 



25 ,2 e-('?ii a ) 2 / 4 
1-S 4 



cos q x L — cosh 



t a La 

" L B 



where Ib = \J h/eB is the magnetic length for a single 
electron. At zero external field ^(qy) takes on the sim- 
plified form of 



A s J m (m,B = 0) 



452 e -( 9|| a) 2 /4 

1-S 4 



• sin 



q x L 

2 



(15) 



Charge distribution of two electrons in a biased 
double quantum dot 



In the case of a singlet-triplet qubit, 1 the DQD is bi- 
ased. The interdot bias is in the regime where the ground 
triplet state remains in the (11) configuration, while the 
ground singlet state is generally a superposition of the 
(11) [denoted as S(ll)] and (02) [denoted as S(02)] sin- 
glets. In S(ll), the two electrons are symmetrically dis- 
tributed across the two dots. In S(02), the two electrons 
are both in the ground orbital state of the lower-energy 
dot. S(ll) and S(02) are tunnel coupled, and the com- 
position of the ground singlet state |S) depends on the 
detuning e between the two singlets in the absence of the 
tunnel coupling. Here e = is defined as the anticrossing 
point of S(ll) and S(02). For negative (positive) e, S(ll) 
[S(02)] has lower energy. 



|S) 
|S(H)> 


= a|S(ll))+/?|S(02)) 
= ^ x -L|t|-|t) 


(16) 
(17) 


|S(02)> 


= \R(1)R(2)) x - It) 


(18) 


|T) 


= Vas x -^=|tl + lt) 


(19) 



(14) 



Here for simplicity we assume a and j3 to be real. Both 
are functions of the interdot detuning e. With the |S) and 
|T) states given and the respective charge distributions 
known, we can calculate the charge distribution differ- 
ence for a biased DQD as a function of a and (3: 



i Biased 



(Hi) 



2 p -(9||a) J /4 



$iTiq x L — a/3 



V2S 



(<?|| a ) 2 / 4 J e «<?*£ _|_ cogn 



q y a La 
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+ 



-(1-S 2 ) 



2S 2 
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"(911 



a f/A | 



q y aLa 
cos q x L - cosh ( -|- 



(20) 
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Similar to the case of a symmetric DQD, at zero magnetic field the expression of is simplified: 

AS 2 



i Biased 



(qii,£ = o) 



-af3^2L e-^^^cos 29 ^ 



VTTs 2 



l-^(l-S 2 ) 



1 - S 4 



(9||°) 2 /4 sm 2^ 



(3 2 +af3 



V2S 



~(?na) 2 /4 sinqxL 



(21) 



The finite interdot bias, which leads to all the addi- 
tional terms in when /3 ^ 0, has some important 
consequences. One distinct feature of Eq. (21) is the 
first term on the right hand side, which does not go 
to zero when q x — > 0. It implies that low- frequency 
phonons are more efficient in causing dephasing for a bi- 
ased DQD compared to a symmetric DQD. In Fig. 1 we 
show A^ lased as a function of q x for various detuning e. 
As discussed above, as e approaches 0, the S(02) compo- 
nent increases in the ground singlet state, and A Blased 
acquires a finite value at q x — 0. On the other hand, 
for e < so that f3 —> 0, the biased DQD system ap- 
proaches the symmetric case, so that A Blased —> A^ vm . 
Furthermore, for the more symmetric DQDs, A$ has a 
peak around q x ~ 1/L, as can be seen from the functional 
form of A^ ym . Another interesting feature of Eq. (21) is 
the last term on the right hand side, which apparently 
does not go to zero when overlap S — » as long as (3 is 
finite. This term is again due to the charge distribution 
difference between (11) and (02) configurations. We will 
explore the consequence of this term at the end of the 
next section. 



D. Two-spin dephasing due to electron-phonon 
interaction with a dissipative bosonic reservoir 




FIG. 1: (Color online) Charge density difference A^ iased as 
a function of phonon wave vector q x and interdot bias e in 
a biased DQD. Notice that as soon as e moves away from 
the S(ll)-S(02) anticrossing point (where e = 0) toward the 
negative bias, A^, has similar characteristics in the form of a 
peak determined by the inter-dot distance L. In this regime 
the two electrons are in the (11) configuration, essentially the 
same as the case of a symmetric DQD. Close to e = 0, the 
two-electron singlet state acquires a (02) component, while 
the peak of A$ shifts toward q x — 0. 



The effective electron-phonon interaction Hamiltonian 
of Eq. (9) is a typical spin-boson Hamiltonian that leads 
to decay in the off-diagonal element of the 2x2 density 
matrix: 68 

PstW = PsT(0)e~ B2(t \ (22) 
where the dephasing factor is positive definite: 68 

^) = -L I d» q |M(q) t (q)|2 - 2 ¥ «** ■ 

'Kit J Uq 2 Kb± 

(23) 

It has long been pointed out that bosonic reservoirs 
with vanishing density of states at low frequencies do 
not cause complete decay of the off-diagonal element of 
a two-level system density matrix, 64 > 69 ~ 77 i n other words 
B 2 (t) of Eq. (23) does not diverge with time, because 
bosonic modes with uj — > determine the long-time be- 
havior for the two-level system. This absence of com- 
plete dephasing can be traced back to the assumptions 
made when the dephasing formula Eq. (23) is derived. 



While it does account for the fact that the bosonic reser- 
voir is in a thermal equilibrium before getting into con- 
tact with the spin, 68 it treats the harmonic modes in 
the bosonic reservoir as completely coherent. However, 
these harmonic modes, in the present case the phonons, 
also belong to an open system, and could lose their co- 
herence to their environments. 78-81 When relaxations of 
the bosonic modes are taken into account, it is expected 
that pure dephasing of the two-level system would even- 
tually become complete. For example, in a spin-boson 
model study of localization, Ref. 82 showed how anhar- 
monicity of the bosonic reservoir, in the specific form of 
two-phonon scattering (which is one out of four terms in 
a transformed phonon interaction Hamiltonian), would 
lead to complete decoherence of the two-level system con- 
sidered. Here we account for phonon relaxation by first 
deriving the phonon Langevin equations 83 that describes 
the effects of phonon-reservoir interactions. 

As we show in Appendix A, the Langevin equation 
for the phonon annihilation operator (in the Heisenberg 
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picture) takes the form 

^Oq(*) = -«WqOq(*) - ySW ~ ^q^ + ^qWe"™"' . 

. (24) 

where w q is the phonon angular frequency, 7q is the 
population relaxation rate of the phonon mode q, 
K q = M\(q)A < f > (q) is the two-electron-phonon interac- 
tion strength from Eq. (9), a z is the Pauli operator in the 
truncated two-electron singlet-triplet space, and F q (i) is 
a noise operator of the reservoir. As shown in Appendix 



A, the specific forms of 7 q and Fq(t) depend on the reser- 
voir and the system-reservoir interaction. However, the 
form of the Langevin equation (24) is quite generic. Com- 
pared to a dissipationless phonon mode, now we have the 
additional second and fourth terms on the right hand 
side of Eq. (24) , representing the dissipation and fluctua- 
tion caused by the coupling to the reservoir. We include 
these two terms and rederive Eq. (23) using the approach 
adopted in Ref. 68. The outline is sketched in Appendix 

B. Now we obtain 



PST(t) = PsT(0)e- B ^)- B ^) > 

V f |M(q)AJq)| 2 f ujI - (7q/2) 2 / \ e-^UalJ* . 1 

27r 3 ft 2 J w q + (7 q /2) 2 I J eF5 a, _ 1 

R 2 m V r 3 |M(q)A (q)| 2 / 7q x 
B * {t) =2^j d q w » + ( 7q /2) 2 iTV = 



At the limit that phonon decay rate 7 q — > 0, Bf(t) — > 
£? 2 (i) while £>f (i) — > 0. For finite 7q , corresponding to 
a dissipative phonon reservoir, we obtain an additional 
exponential decay of the off-diagonal density matrix ele- 
ment in Eq. (26) compared to the non-dissipative reser- 
voir result of Eq. (23). The rate of this exponential decay 
Tst is proportional to the phonon decay rate 7 q inte- 
grated over the phonon modes. Notice that T$t docs 
not explicitly contain the thermal factor coth that 
describes the thermal occupation of the phonon modes. 
This is because B 2 (t) comes from phonons decaying into 
their reservoirs, when phonons themselves are regarded 
as coherent bosons, while temperature information of the 
reservoirs for the phonons is contained in the decay rate 
7q and the noise operator F q (£). 



III. RESULTS 



The main questions we would like to answer in this 
work are as follows: Is electron-phonon interaction an 
important decoherence channel for spin qubits in semi- 
conductor quantum dots? Under what condition is it 
important? How do different substrate materials (GaAs 
and Si) compare to each other? And how do different 
qubit architectures compare with each other? Below we 
show our results that provide the answers. 



A. Symmetric double dot 

Let us first examine the dynamical behaviors of the 
dephasing factors B 2 (t) and dephasing rate Tst due to 
electron-phonon interaction when the double quantum 
dot is unbiased. 

In Fig. 2 we show the typical behavior of the dephas- 
ing factor B 2 (t) in the absence of phonon decay for var- 
ious types of electron-phonon interactions in GaAs and 
Si. There are two interesting features shared by all the 
curves for acoustic phonons in Fig. 2. At very short times 
ps), the increase of B 2 (t) is quadratic, which orig- 
inates from Taylor expansion of the sin 2 Wqi/2 factor in 
the integrand at the small-i limit. At long times, all the 
curves saturate, which means that dephasing does not 
increase with time anymore, so that it corresponds more 
to a finite loss of contrast than the conventional complete 
decay of off-diagonal density matrix elements. The tran- 
sition between the quadratic increase and the saturation 
happens between 1 and 10 ps for double dots with a dot 
separation of 40 nm and a single-dot radius of 20 nm be- 
cause this time is essentially determined by the interdot 
distance divided by the speed of sound (5~8x 10 3 m/s 
in Si and 3 — 5 x 10 3 m/s in GaAs): 40 nm / c — 10 
ps. The saturation time for Si is shorter because Si has a 
larger speed of sound. The dephasing factor due to DP 
interaction with TA phonons in Si is two orders of mag- 
nitude smaller than that due to LA phonons, and is not 
plotted in Fig. 2. 

Mathematically, the long-time saturation can be un- 
derstood by writing 2sin 2 w q i/2 as 1 — coswqi = 1 — 



cos(cqt). Since the acoustic phonon spectrum is contin- 
uous and nonsingular, the cosine term leads to a van- 
ishing contribution to the integral at large times, which 
leaves the dephasing factor determined by a constant in- 
tegral that is independent of time. Physically, this sat- 
uration is due to the fact that long-time dephasing is 
determined by the low-frequency part of the spectrum of 
the bosonic reservoir, while the phonon density of states 
vanishes quadratically at low frequency. In other words, 
non-dissipative acoustic phonons simply form an ineffi- 
cient dephasing reservoir as compared to other charge 
fluctuation reservoirs such as fluctuating charge traps, 
which have a 1/f spectral density. This incomplete de- 
phasing has been observed theoretically in a variety of 
calculations related to phonons, in the studies of general 
spin-boson decoherence behaviors, charge and spin coher- 
ence of single electrons, and exciton coherence. 64,69-77 

For electron interaction with optical phonons in GaAs, 
the dephasing factor B 2 (t) takes on a particularly simple 
form because the optical phonon dispersion at the zone 
center is flat. Take uj„ rs u)LOi we obtain 



Bloit) 



V f , 3 jM(q)A (q)| 2 ^ 2 . q , ^ 



ir 3 H 2 
2e 2 



efqi 
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TT 2 hlJLO \eoo 

d a q l^(q)l 2 



■ sin 



k B T 



1 \ , i fiWLO . 2 u LOt 

— coth sin 

eoj k B T 2 



b 2 LO sm 



i 

2 ^Lot 



(27) 



which is a sinusoidal function of time. For a GaAs double 
dot with a single-dot radius of a = 20 nm, and L/a in 
the range of 1 and 2, the coefficient b 2 LO for the sinusoidal 
function ranges between 10 -4 and 10~ 9 . 

As we discussed in the previous section, the abso- 
lute value of the saturated dephasing, as long as it is 
small (<C 1), is not an important parameter by itself 
because dephasing will eventually become complete due 
to phonon relaxation. However, the relative magnitudes 
of the saturated dephasing shown in Fig. 2 do give a 
qualitative sense of the relative importance of various 
types of electron-acoustic-phonon interactions. Specifi- 
cally, in GaAs PE coupling to TA phonons produces the 
strongest dephasing effect, while in Si DP coupling to 
LA phonons is the most important. In addition, as indi- 
cated in Eq.(23), B 2 {t) does have a strong temperature 
dependence as well. At higher temperatures more acous- 
tic phonon modes contribute to dephasing, so that B 2 (t) 
can eventually become an O(l) quantity and dephasing 
can be considered complete. In Fig. 3 we plot the tem- 
perature dependence of the saturated B 2 (t) for GaAs and 
Si quantum dots. At temperatures above 1 K dephasing 
increases with temperature almost linearly. On the other 
hand, optical phonon induced dephasing does not have a 
pronounced temperature dependence even at T = 10 K 
because huj^o ~ 36 meV is much larger than k B T at low 
temperatures. 
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FIG. 2: (Color online) Two-spin dephasing in a symmetric 
double quantum dot induced by a non-dissipative phonon 
reservoir. All the curves are for double quantum dots with 
an interdot separation of 40 nm and single dot orbital radius 
of 20 nm. More specifically, the black solid line is for PE 
interaction with TA phonons in GaAs, the red dotted line is 
for PE coupling to LA phonons in GaAs, the green dashed 
line is for DP coupling to LA phonons in GaAs; the blue dot- 
dashed curve is for DP coupling to LA phonons in Si; and 
the maroon dot-dashed-dashed horizontal line represents the 
dephasing magnitude for polar interaction with LO phonons 
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FIG. 3: (Color online) Phonon induced two-spin dephasing 
rate as a function of acoustic phonon temperature in a sym- 
metric double dot for both GaAs and Si. The single dot wave 
function radius for all the data is 20 nm. 



When phonon decay is included, the most important 
effect is the added exponential dephasing e _rsT *, upon 
which we will focus in the rest of this study. To calcu- 
late Tst, we need to first identify the q-dependence of 
the phonon relaxation rate 7 q . Qualitatively, lower en- 
ergy (lower frequency w q ) acoustic phonons have to decay 
slower (smaller 7 q ). For example, when phonon decay is 
due to anharmonicity, or more specifically the third order 
process of one phonon splitting into two, the phonon de- 
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cay rate could vary as q n with n between 1 and 4 depend- 
ing on the lattice symmetry and phonon branches. 78 ' 79 In 
the case of a phonon cavity, it is not clear how the Q- 
factor would vary with the phonon wave vector, although 
with 7 q oc -D(wq) [where -D(w q ) is the phonon reservoir 
density of states] one could expect 7 q oc q 2 . In the follow- 
ing, we calculate the acoustic-phonon-induced dephasing 
rate assuming that 7 q = 709", with n taking the value of 
2 or 3. Taking Q = 10 3 for a TA phonon with an energy 
of 0.1 meV, we obtain 70 = 10 8 1/s. For LA phonons, 
which have higher energies than TA phonons with the 
same q, we take 70 = 10 9 1/s. This is an arbitrary choice 
that is used to reflect the fact that LA phonons gener- 
ally have shorter lifetimes than TA phonons. 79 For LO 
phonons we assume a constant relaxation time of 10 ps 
for all modes. The rationale here is that LO phonons 
have a flat dispersion, so that they should have a near 
constant relaxation rate near the Brillouin zone center. 
The LO phonons also have a very short lifetime because 
of their large energy. Zone center LO phonons have been 
measured to have a lifetime of 7 ps. 84,85 

In Fig. 4, we plot the phonon-induced two-spin de- 
phasing rate Tst in symmetric DQDs in both GaAs and 
Si as functions of the half interdot distance L. The ra- 
dius of the single dot electron wave function is 20 nm for 
this figure and all the following figures. The strong de- 
pendence on L for all data sets originates from the fact 
that the charge distribution difference between the two- 
electron singlet and triplet states in a symmetric DQD 
is directly dependent on interdot wave function overlap: 
Tst oc [4S' 2 /(1 — S* 2 )] 2 , so that the smaller the over- 
lap, the smaller the difference in charge distribution, and 
the smaller the phonon-induced dephasing. Based on the 
data given in Fig. 4, phonon-induced dephasing is not an 
important decoherence mechanism when L/a > 2 in a 
symmetric DQD. 

An important feature of Fig. 4 is that DP coupling 
to LA phonons is the most important dephasing chan- 
nel for GaAs, and produces about the same magnitude 
of dephasing in Si. In GaAs, dephasing due to DP cou- 
pling is about one order of magnitude larger than that 
by PO coupling to LO phonons, and almost two orders 
of magnitude larger than that due to PE coupling to 
both LA and TA phonons. This fact is somewhat sur- 
prising because in Fig. 2 it is clear that PE coupling to 
TA phonons is by far the most important decoherence 
channel. However, notice that in the present calcula- 
tion of Pst the acoustic phonon decay increases rapidly 
as phonon energy increases, so that the contributions 
from higher-energy phonons are much more important 
in the calculation of Tst than in B 2 (t). This tilt to- 
ward higher-energy phonons strongly favors DP coupling 
over PE couplings because of the factor of q difference in 
the electron-phonon coupling matrix element. The simi- 
lar values for dephasing for GaAs and Si within the DP 
mechanism is more of a coincidence: they have similar 
values in DP coupling strength, mass density, and speed 
of sound, and we chose the same 70 for both materials, 



although Si has the extra contribution from the shear 
DP constant S u . Based on the results presented in this 
figure, phonon-induced dephasing is an essentially equiv- 
alent decoherence mechanism for Si and GaAs. 

Another interesting aspects of Fig. 4 is that LO 
phonons turn out to be a strong source of dephasing for 
the two-spin states in GaAs, even though they have very 
high energy in GaAs (~ 36 meV). This somewhat surpris- 
ing result originates from the facts that LO phonons have 
a diverging density of states at the zone center (as com- 
pared to the vanishing density of states for the acoustic 
phonons) and a very fast relaxation rate (experimentally 
measured at 7 ps 84 ), and that GaAs has a reasonably 
strong polar interaction strength. In Si, the conduction 
electrons do not interact with optical phonons, therefore 
this decoherence channel is completely removed. 

For DP and PE interactions, the different in- 
dependence of the phonon relaxation rate 7 q leads to 
quite different results in the two-spin dephasing rate Tst- 
With DP interaction, Tst is not very sensitive to the ex- 
ponent n, and increasing n in 7 q ~ q n leads to a slight 
decrease of Tst- On the other hand, for PE interaction, 
increasing n leads to an approximately three-fold increase 
of Tst- The change of the exponent n leads to a shift 
of the dominant q region that contributes to dephasing. 
For PE interaction, increasing n from 2 to 3 moves the 
dominant contribution to larger q phonons, which have 
a larger density of states, leading to an increase in Tst- 
For DP interaction, the dominant contribution already 
comes from the qa ~ 1 region, where changing n does 
not have much of an effect. 

In Fig. 5, we plot the two-spin merit figure M. as a 
function of the interdot distance for double dots in GaAs. 
Here the merit figure is defined as the ratio between a 
typical exchange gate time given by h/J (J is the ex- 
change splitting) and the two-spin decay time given by 
1/rsT: M. = J /fiT 'st- The exchange splitting J is cal- 
culated within the Hcitler-London model with a quar- 
tic confinement potential 86 . The increase of the merit 
figure at larger inter-dot distance reflects the fact that 
the exchange splitting and the phonon-induced dephas- 
ing have a different dependence on the interdot overlap 
integral S: J <~ S 2 , while Tst ~ S 4 . The results shown 
in this figure indicate that for a two-dot exchange gate 
to operate with a low error rate, a slower operation with 
a smaller interdot overlap is preferable with regard to 
phonon-induced dephasing, and fault-tolerant two-qubit 
operations should be achievable for pretty strongly cou- 
pled dots, with L/a > 1.5. We do not have any data for Si 
DQDs in this figure. Calculating exchange interaction in 
a Si double dot requires much more sophisticated quan- 
tum chemical approaches than a simple Heitler-London 
approximation 51,87 ' 88 because in Si the interaction effect 
is stronger compared to GaAs (larger effective mass and 
smaller dielectric constant), so that the Heitler-London 
approximation does not adequately account for the two- 
electron correlation. For the current evaluation, it is 
sufficient to point out that Fig. 4 above indicates that 
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FIG. 4: (Color online) Phonon-induced two-spin dephasing 
rate as a function of the interdot separation for a GaAs and 
a Si symmetric double quantum dot. The horizontal line is 
drawn at a dephasing time of 1 /is, approximately the de- 
coherence times measured in Refs. 1,2. The legends for the 
data sets have the following format: type of materials (GaAs 
or Si), type of interaction (DP, PE, or PO), type of phonons 
involved (LA, TA, or LO), and the g-dependence of 7 q (q 2 , 
q 3 , or constant 7_lo). 



phonon-induced dephasing is about the same order of 
magnitude in Si as in GaAs, while exchange coupling 
should only be somewhat smaller than in GaAs. There- 
fore overall the merit figure should remain about the 
same when moving from GaAs to Si. 



B. Biased double dot 

For a biased DQD, the main question is whether the 
admixture from the S(02) singlet state and the result- 
ing dipole coupling would lead to significantly increased 
dephasing. Interestingly, the bias not only directly af- 
fects the value of A^(qii), but also its functional form. 
In Fig. 1 we have shown how A<j>(q x , q y = 0) depends on 
q x for various interdot biases e. It is clear from that fig- 
ure that the peak of shifts toward q x = 0, while the 
peak height decreases, as the interdot bias shifts from 
the (11) toward the (02) regime. Furthermore, Eq. (21) 
indicates that as soon as f3 ^ 0, there is a mixing of S(ll) 
and S(02) states, so that A^ tased acquires a nonvanishing 
component [1st term on the right hand side of Eq. (21)] 
as q x — > 0, leading to an increase in the phonon-induced 
dephasing. 

Now we can calculate the two-spin dephasing rate Tst 
for any voltage bias between the dots. Figure 6 shows 
Fst as a function of dimensionless interdot bias e. As e 
becomes increasingly negative, the biased DQD states ap- 
proach those of a symmetric DQD, and T$t approaches 
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FIG. 5: (Color online) Merit figure based on phonon-induced 
dephasing of two-spin states in a symmetric GaAs double dot 
as a function of half interdot distance. We draw a line at 10 4 
as the nominal threshold for fault tolerant quantum computa- 
tion. Therefore the double dot (with single-dot wave function 
radius at 20 nm) should be kept apart further than 60 nm. 
The legend format is similar to in Fig. 4 (without the first 
item for materials as all data here are for GaAs): type of 
coupling, type of phonon, and q-dependence of 7 q . 

the value given in Fig. 4. On the other hand, as e in- 
creases toward positive bias, the ground singlet state has 
a larger S(02) component, and A$ a larger dipolar contri- 
bution, so that F"st increases. At e = 0, Tgx is dominated 
by the dipolar contribution from the S(ll)-S(02) mixing, 
and is about ten times larger than in a symmetric DQD, 
where A,p is determined by a quadrupolar charge distri- 
bution difference between S(ll) singlet and T(ll) triplet 
states. 

As indicated in Fig. 4, for symmetric DQDs phonon- 
induced decoherence becomes much less important at 
larger L because of the overlap factor S in the charge 
difference A$. In the case of a biased DQD, when the 
DQD has a vanishing overlap, 

Af ased {B = 0, S = 0) = -z/3 2 e- (9 n a)2/4 sin q x L , (28) 

which does not seem to depend on the interdot over- 
lap. This term leads to dephasing between T(ll) and 
S(02) states, which clearly have different charge distri- 
butions. However, if overlap S vanishes because the in- 
terdot distance L increases, phonon-induced dephasing 
will not saturate to a constant, as Eq. (28) seems to in- 
dicate, because /3 depends on L as well. Recall that the 
tunnel coupling t of the S(ll) and S(02) singlet states 
is t = (S(11)|JT|S(02)) oc S = e- L2 / a \ When L in- 

/ r 2 r 2 \ / 2 

creases, t decreases as t = to e~ ( " « w ° for two cou- 
pled parabolic dots, where to is the tunnel coupling at Lq. 
For a fixed interdot bias e [recall that e = corresponds 
to the S(ll)-S(02) crossing point, where /3 2 = 0.5], the 
DQD moves farther from the anticrossing point as t gets 
smaller, leading to a decreasing (3 [which is the weight of 
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FIG. 6: (Color online) Phonon-induced two-spin dephasing 
rate as a function of the interdot bias e for a GaAs dou- 
ble quantum dot with piezoelectric coupling to transverse 
phonons. The double dot are separated by 40 nm and the 
single-dot wave function radius is 20 nm, so that we are at the 
strong coupling limit. At large negative bias the relaxation 
rate approaches a value of about 10~ 6 1/ps, or T$ of about 1 
/is. As the bias increases toward the S(ll)-S(02) anticrossing, 
the dephasing rate increases so that at e = 0, ~ 100ns. 
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FIG. 7: (Color online) Phonon-induced two-spin dephasing 
rate as a function of the interdot distance L for a biased GaAs 
double quantum dot with DP coupling to LA phonons, PO 
coupling to LO phonons, and PE coupling to TA phonons. 
The format of the legend is the type of coupling, the type of 
phonon, and the g-dependence of the phonon relaxation rate 
7 q . The single-dot wave function radius is 20 nm, and the 
interdot bias is e — — 1 at L = 20 nm. 



the higher-energy singlet; for negative e, it is the weight 
of the S(02) state]. In Fig. 7 we plot the dephasing rate 
for a biased DQD as a function of the interdot separation. 
The figure shows the same rapid decrease of dephasing 
for all types of phonons as L increases, similar to the sit- 
uation in symmetric DQDs. Indeed, putting ft oc S into 
Eq. (21), it is clear that A b j ased oc S 2 , which is the same 
overlap-dependence as in the case of a symmetric DQD. 

The results in Figs. 6 and 7 show that in a biased 
DQD phonon-induced dephasing approaches the sym- 
metric DQD limit at large negative bias, and increases 
monotonically as interdot bias e increases. However, one 
can always reduce this dephasing by increasing the in- 
terdot distance and reducing the wave-function overlap. 
Furthermore, for larger negative biases [deeper into the 
(11) regime, with smaller exchange splitting J = i/2e], 
/3 2 is smaller. When |e| > 1, /3 2 ~ (1/e) 2 while J oc 1/e. 
Therefore, there should exist a regime where the dephas- 
ing rate Tst is much smaller than exchange splitting, 
so that fault-tolerant exchange gates can be performed. 
For example, if we choose L/a = 2 with a = 20 nm, 
l/r sx ~ 100 fxs even at e = —1, according to Fig. 7. 
At such an interdot separation, \t\ ~ 10/ieV, so that 
J ~ 1 //eV for e — —10, with a gate time in the order of 
a nanosecond, leading to a merit figure of ~ 10 5 . 

IV. DISCUSSIONS AND CONCLUSIONS 

Based on our results presented in this study, phonon 
induced two-spin dephasing in both symmetric DQD's 
and biased DQD's can be strongly suppressed by reduc- 



ing the double-dot tunnel coupling. The strong overlap 
dependence of the dephasing rate dictates that phonon- 
induced dephasing is only important when the double 
dot is tightly coupled. Dephasing for a biased double 
dot does increase with bias because of the admixture of 
the S(02) state in the singlet ground state, which intro- 
duces electric dipole coupling into phonon-induced deco- 
herence. Therefore, phonon-induced two-spin dephasing 
is generally stronger in biased DQD's, such as in the case 
of a singlet-triplet qubit. 

Phonon-induced two-spin dephasing studied in this pa- 
per is related to the different dressing that singlet and 
triplet electronic states experience through interaction 
with the phonons. When phonons themselves decohere, 
this spin dephasing channel leads to true complete de- 
coherence. On the other hand, ensemble average over 
phonon modes (while each evolve coherently) only leads 
to a finite degree of dephasing. A legitimate question here 
is whether this part of the dephasing (due to phonon pop- 
ulation average) would disappear if we consider dressed 
electron spin states, especially considering that this fi- 
nite dephasing generally saturates in the order of 10 ps, 
much faster than the electron spin initialization and ma- 
nipulation processes in quantum dots. Mathematically, 
the answer to this question may very well be "yes" , as 
long as one can identify the energies of the dressed states 
precisely. But to answer this question with confidence, 
one needs to clarify how the energies of the spin states 
are measured, and how electron-phonon interaction may 
be incorporated in the description of measurement. In 
the current generation devices, spin detection is achieved 
through charge sensing in the spin-blockade regime, 89 
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which is insensitive to phonons, so that the pho non- 
induced dephasing due to ensemble averaging cannot be 
removed. Ultimately, though, this question is moot be- 
cause phonons do relax and are not coherent forever. 

The phonon-induced pure dephasing mechanism we 
consider here originates from the charge distribution dif- 
ference between states that have different spatial sym- 
metry, and involves no real or virtual phonon emission 
or absorption. It is different from another mechanism 
of phonon-induced dephasing studied in Ref. 54, which 
is based on different level distribution of electron singlet 
and triplet states and involves virtual emission and ab- 
sorption of phonons. 

In conclusion, we have studied phonon-induced pure 
dephasing between two-electron singlet and triplet spin 
states in a semiconductor double quantum dot. We find 
that this pure dephasing is important for tightly coupled 
double dots, but is strongly suppressed when the dou- 
ble dot separation increases, so that at relatively large 
dot separations (L/a > 2) fault-tolerant exchange gates 
can be realized. A biased double dot has stronger de- 
phasing compared to a symmetric double dot with the 
same dot parameters due to the mixing of (11) and (02) 
singlet states and the resulting finite electric-dipole cou- 
pling. We have quantified two-spin dephasing in both 
GaAs and Si double dots, finding that deformation po- 
tential coupling to LA phonons is the most important 
dephasing mechanism in both materials, and produces 
about the same magnitude dephasing in both materials. 
We also find that the LO phonon makes a non-negligible 
contribution to dephasing in GaAs because of the very 
fast optical phonon relaxation. Overall, phonon-induced 
two-spin dephasing is an equivalent decoherence mech- 
anism for Si and GaAs, is stronger in a biased double 
dot than in a symmetric double dot, and can be sup- 
pressed by reducing the interdot overlap of the electron 
wave functions. 
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from useful discussions with Peter Yu, Luming Duan, 
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Appendix A: Phonon relaxation 

We can model phonon relaxation by assuming that 
each phonon mode couples to a continuum of bosonic 
modes. This is identical to the description of a cavity 
photon mode coupling to a continuum. 83 Such a model is 
grounded in the development of phonon cavities in semi- 
conductor heterostructures, 80 ' 81 but more importantly, 
it presents a clear physical picture about phonon relax- 



ation. As we discuss below, the exact form of phonon- 
reservoir interaction does not change the general features 
of phonon relaxation. The Hamiltonian for a phonon 
mode and its reservoirs takes the form 

H = H s + H r + H int 
H s = hjj^a^ay, 

q 

H int = [-9( k ' °l) a k & q + 5*(k,q)a k & q , 

q 

where is the phonon annihilation operator, 6 q is the 
annihilation operator of the bosonic modes in the reser- 
voir, and g(k, q) is the coupling strength between the 
phonon mode and the reservoir modes. With this inter- 
action with the reservoir, the Langevin equation for the 
phonon annihilation operator (in the Heisenberg picture) 
takes the form (as is obtained in the discussion of cavity 
photon decay in any number of quantum optics books, 
e.g., Ref. 83, Chap. 14, Sec. 14.3) 

j t a k (t) = -iw k a k (t)- ^a k (t)+F k (i)e- iUkt (Al) 

F*(t) = -^ 5 (k,q)6 q (io)e-^^- Wk)t (A2) 
q 

7k = 27rD(w k ) |<?K)| 2 , (A3) 

where 6 q = 6 q e* Wq * is the slowly varying amplitude of the 
reservoir bosonic operator. Notice that the noise opera- 
tor Fk(£) here is assumed to be independent of the initial 
time to. In the definition of the decay rate 7 k the sum 
over reservoir mode q has been replaced by an integra- 
tion over energy, with D(uj) being the density of state for 
the reservoir, and g(k, q) is assumed to be smooth in any 
narrow energy range so that it can be replaced by g(u>). 

In a realistic system, phonon relaxation often orig- 
inates from phonon anharmonicity. The lowest-order 
anharmonic interaction has each phonon mode coupled 
to pairs of other modes in the phonon reservoir (three- 
phonon processes) . We can choose any particular phonon 
mode as the system and consider its dynamics coupling 
to the rest of the phonon modes, 

H = Hq + Hint 

Ho = fiw q a q a q 

Hint = [.9(q>q') a q a q' a q-q' +.9*(q,q') a q a q' a q- q ' 

q,q' 

For simplicity and clarity, we have limited ourselves to 
normal processes for phonon relaxation. 66 Following the 
same standard procedure as above, 83 we first obtain the 
Heisenberg equations for the operator a q , then formally 
solve its equation of motion, and finally put the results 
back into the equation of motion for a q . The only sig- 
nificant difference here is that each phonon mode plays 
both the role of a higher energy phonon relaxing to two 
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lower energy ones, and the role of a lower energy mode to the lowest order contributions, we obtain 
which a higher-energy phonon can decay. Keeping only 



-a q (i) = -iu^t) - fa^t) + F^t)e-^ (A4) 

*q(*) = -tY,9*(^^(to)~a^(to)e-^' + ^'-^ - i £ g(q[ , q^^-q^e"^-^-.,-^)' (A5) 
q q 

7q = ^2 2 \9(^q!)\ (1 +«q' +«q-q') (5 ( W q ~ W q' - W q-q' ) + X] 2 l.9(q' > q) I 2 («q'-q - n q /)5(w q / - W q - Wq'-q) ,(A6) 

q' q' 

I 



where a q = a q e lWc| * is again the slowly varying amplitude 
of the phonon operator. In obtaining the Langevin equa- 
tion (A4) for the phonon operator a q we have assumed 
that the population of the relevant reservoir phonon 
modes is not strongly perturbed by the phonon inter- 
action, which allows us to factor out the phonon oper- 
ators from the time integrals and replace number op- 
erator a q a q by a constant thermal phonon population 
number n q . Notice that the expression for 7 q contains 
two parts. The first part is due to two-phonon emission 
and is proportional to 1 + n q / + n q _ q / , including both 
spontaneous and stimulated emissions. The expression 
is consistent with the population decay rate obtained 
using Fermi's Golden Rule, which is also proportional 
to 1 + riqi + ra q _ q '. 66 The second part of 7 q is due to 
the absorption process where the q-mode phonon and a 
(q' — q)-modc phonon are converted into a higher-energy 
q'-mode phonon. This contribution is always positive at 
thermal equilibrium since n q / < n q /_ q . Another inter- 
esting feature is that while in the case of cavity phonon 
decay 7 q is only determined by the reservoir density of 
states and coupling strength [see Eq. (A3)], here it is 
also influenced by the thermal occupation of the reser- 
voir phonon modes. 

An important result above is that Eqs. (A4) and (Al) 
are of exactly the same form. This is not really surprising: 
both coupling to a bosonic reservoir through a cavity mir- 
ror and anharmonic interaction with other bosonic modes 
lead to relaxation, as has been clearly demonstrated in 
numerous experiments and clarified in many theoretical 



studies. 66 ' 83 We can therefore use Eq. (24) to describe 
dissipative phonon dynamics, with appropriately chosen 
phonon relaxation rate 7 q and noise operator Fq(t). 

Appendix B: Phonon induced spin dephasing 

The effective spin-boson Hamiltonian for the two- 
electron-phonon interaction is given by Eq. (9): 

Heff = ^M A (q)^(q)a z (a q , A + aL qjA ) . (Bl) 

This is a pure dephasing Hamiltonian: the interaction 
with phonons does not change the two-electron eigen- 
states and there is no relaxation due to phonons. We 
derive the phonon induced dephasing following the ap- 
proach presented in Ref. 68. The first step is to obtain 
the complete phonon Langevin equation (the Heisenberg 
equation in Ref. 68 because phonon relaxation is not con- 
sidered there), which contains both the terms as derived 
in Appendix A and a term from the electron-phonon 
interaction Hamiltonian above. The phonon Langevin 
equation now takes the form 

(B2) 

where K q = M\(q)A ( f ) (q) for any particular phonon 
branch. This differential equation for the phonon op- 
erators can be formally solved. The solution is 



Oq(t) = a q (0)e- ( ^ +iw ^* + ■. 



ioj. 



+ la 

^ 2 



+ 



/' 

J a 



dt'Fq(t')e 



(B3) 



This solution is formal because the noise operator F q (t) 
contains a sum over all the phonon operators themselves. 
However, when we consider a phonon reservoir that is 
not driven far away from its equilibrium, which is the 



case studied in this work, the noise operator would then 
essentially be a mean field average over all the phonon 
modes and can be treated as independent from individual 
phonon mode properties. Therefore, to a good approxi- 
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mation, Eq. (B4) gives the phonon evolution in the pres- 
ence of electron-phonon and phonon-phonon or general 
phonon-reservoir interactions. 

The time evolution of the phonon operators allows 
us to solve for the total density matrix of the electron- 
phonon system, because the electron-phonon interaction 
is diagonal (i.e. pure dephasing spin-boson interaction), 
as discussed in Ref. 68. The phonon operators can then 
be traced out using a coherent-state representation for 
the phonon modes, assuming an initial thermal equilib- 
rium distribution. 68 One subtle point in the derivation 
of the total density operator is that there is a sign dif- 



ference between the von Neumann equation for density 
operators and the Heiscnbcrg equation for regular oper- 
ators. While one can expand a density operator in terms 
of regular operators, p(Q) — J2 n a nA n (0), the time evolu- 
tion is reversed: p(t) = ^2 n a n A n (—t). Thus we need to 
calculate a q (— t) before calculating the density operator. 
It is important to note here that for a q (— t), the phonon 
decay term e~ 7qt / 2 remains a decay term, because 7 q 
comes from a second-order perturbation calculation and 
is not affected by the change of time direction: 



= a q (0)e" 



-(2a_; Wq)t 



_(-^._i Wq )t _ 
-«"q + ¥ 



dt'F^{-t')e 



-(^- iWq )(t-t') 



(B4) 



The key difference between the present calculation and 
the one in Ref. 68 is that here the solution to the phonon 
operator contains the additional terms for phonon decay 
and noises. The decay term is fully integrated in our 
calculation and leads to the exponential decay of the off- 
diagonal matrix element of the two-spin density matrix. 
The noise term, as we discussed before, represents an 



overall mean field noise that is independent of individual 
phonon modes or the electrons. With the assumption 
that the noise term commutes with the system density 
operator (i.e. assuming the noise is classical), its effect 
would be to induce a phase shift in the evolution of the 
reduced electron density operator but not decay. 
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